Harmonic oscillator model for current- and field-driven magnetic vortices 
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In experiments the distinction between spin-torque and Oersted-field driven magnetization dynamics is still 
an open problem. Here, the gyroscopic motion of current- and field-driven magnetic vortices in small thin- 
film elements is investigated by analytical calculations and by numerical simulations. It is found that for small 
harmonic excitations the vortex core performs an elliptical rotation around its equilibrium position. The global 
phase of the rotation and the ratio between the semi-axes are determined by the frequency and the amplitude of 
the Oersted field and the spin torque. 
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Recently it has been found that a spin-polarized current 
flowing through a magnetic sample interacts with the mag- 
netization and exerts a torque on the local magnetization.''- 
A promising system for the investigation of the spin-torque 
effect is a vortex in a micro- or nanostructured magnetic thin- 
film element. Vortices are formed when the in-plane mag- 
netization curls around a center region. In this few nanome- 
ter large center region-, called the vortex core, the magneti- 
zation turns out-of-plane to minimize the exchange energyii 
It is known that these vortices precess around their equilib- 
rium position when excited by magnetic field pulses^ and it 
was predicted that spin-polarized electric currents can do the 
samcii The spacial restriction of the vortex core as well as 
its periodic motion around its ground state yield an especially 
accessible system for space- and time-resolved measurements 
with scanning probe and time-integrative techniques such as 
soft X-ray microscopy or X-ray photoemission electron mi- 
croscopy .^'^'^'^''° Magnetic vortices also occur in vortex do- 
main walls. The motion of such walls has recently been inves- 
tigated intensivelyJii^ Understanding the dynamics of con- 
fined vortices can give deeper insight in the mechanism of 
vortex-wall motion. An in-plane Oersted field accompany- 
ing the current flow also influences the motion of the vortex 
core. For the interpretation of experimental data it is crucial 
to distinguish between the influence of the spin torque and of 
the Oersted fieldiii 
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FIG. 1 : (a) Scheme of the magnetization in a square magnetic thin- 
film element with a vortex that is deflected to the right, (b) Magne- 
tization of a vortex in its static ground state. The height denotes the 
z-component while the gray scale corresponds to the direction of the 
in-plane magnetization. 



In this paper we investigate the current- and field-driven gy- 
roscopic motion of magnetic vortices in square thin-film ele- 
ments of size / and thickness t as shown in Fig.[T]and present 
a method to distinguish between spin torque and Oersted field 
driven magnetization dynamics. In the presence of a spin- 
polarized current the time evolution of the magnetization is 
given by the extended Landau-Lifshitz-Gilbert equation 
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with the coupling constant bj = P hb / [eMs{l + between 
the current and the magnetization where P is the spin polar- 
ization, Ms the saturation magnetization, and ^ the degree of 
non-adiabaticity.'^ If the vortex keeps its static structure, its 
motion with the velocity v can be described using the Thiele 
equationji^ This equation was expanded by Nakatani et alJ^ 
to include the action of a spin-polarized current flowing in the 
sample. 



F + G X {v + b-jj) + D{av + ^bjj) = 0. 
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Denoting the out-of-plane angle of the magnetization with 9 
and the angle of the in-plane magnetization with (f>, the force 
due to the external and the stray field is 



F = -^io 

The gyrovector 



dV 



{Hs, ■ M). (3) 



G = - 



f dV sm{e){V0 X V(/.) 
7 J 

2T:MsfJ.otp - 



(4) 



7 



Go 



indicates the axis of precession and points out-of-plane. The 
dissipation tensor is given by 
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It is diagonal with 
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The constant a is the lower bound of the integration. It is in 
the order of magnitude of the radius of the vortex core. -^i'-^''^-'^ 
A polarization p of +1 (—1) denotes that the magnetization 
in the vortex core is parallel (antiparallel) to the z-axis. The 
velocity of the vortex core is in-plane and hence perpendicular 
to the gyrovector. Thus Eq. (|2]) can be rewritten as 

GxF-Gl{v + bjj) + DoGx {av + £,bjj) = 0. (7) 

By calculating G x v from Eq. (|7]) and inserting the result in 
Eq. (|2]i we can derive the velocity 
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of the vortex core. As for any square-symmetric confining 
potential, the stray-field energy for small deflections can be 
modeled as a paraboUc potential 

1 



with the coordinates X and Y of the vortex core (see Fig.[T^). 

In the following a spacially homogeneous current in x- 
direction is investigated. Due to possible inhomogeneities in 
real samples the current flow may vary in the out-of -plane di- 
rection. This results in an in-plane Oersted field which is per- 
pendicular to the direction of the current flow. In the following 
this Oersted field is accounted for by a homogeneous magnetic 
field in y-direction. Both driving forces may depend on time. 
To estimate the Zeeman energy due to the Oersted field H, 
the magnetization pattern is divided into four triangles (see 
Fig.[T^). Assuming that the magnetization is uniform in each 
of these triangles the total Zeeman energy is given by 
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with the chirality c of the vortex. A chirality of +1 ( — 1) de- 
notes a counterclockwise (clockwise) curling of the magne- 
tization around the vortex core. We will see that this simple 
approximation describes the field-induced vortex motion suf- 
ficiently well. In this case the force is given by 

F = —^{Eg+Ez) = —po^IsHltcex—rnuj^Xex—muj'^Yey. 

(11) 

Inserting Eq. ( fTTT i in Eq. ([8]) yields the equation of motion for 
the vortex. In the absence of current and field the excited vor- 
tex performs an exponentially damped spiral rotation around 
its equilibrium position with its free frequency 
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and the damping constant 
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From Eqs. ( fT2] ) and ( fTSl ) one easily obtains that 

Doa = . 



(14) 



For thin-film systems {t/l < 0.1) the resonance frequency of 
a vortex is proportional to the inverse lateral dimension 1 /L— 
Here, we obtain from Eq. ( fT4] i that the damping constant F 
also has a characteristic length dependence, F oc \n{l/a)/l. 
Substituting Doa using Eq. ( fT4l i the equation of motion of the 
vortex can be written as 



-r 

^puj 



—puj 

-r 



hr2 



Go 



Go 



(15) 



— a 1 ■ 
0,1 



In the following we assume harmonic excitations, i.e., the 
magnetic field and the electrical current are of the form 
H{t) = i/oe'"* and = joe'^^*. The magnetic (Oersted) 
field and the electrical current are in phase. Assuming that 
the squared Gilbert damping is small (a^ <C 1), the damp- 
ing constant of the vortex is small compared to its frequency 



(9) (F^ < u?). Then Eq. ^ has the solution 
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with H = 'yHolc/{2TT) and j ~ bj jo- The first two terms 
with prefactors A and B are exponentially damped and depend 
on the starting configuration. Independent of the source of 
excitation, i.e., field or current, the sense of rotation of the 
vortex is given by its polarization, i.e., p = +1 (p = — 1) 
denotes a counterclockwise (clockwise) rotation of the vortex 
core. Changing the sign of the chirality has the same effect 
as turning the magnetic field by 180°. Similar to the motion 
of magnetic domain walls in thin nanowires2i the vortex is 
driven by the current and the magnetic field as well as by their 
time derivatives. 

At resonance the amplitude of the vortex core displacement 
in X- and i/-direction is the same and the vortex performs a cir- 
cular rotation. A vortex which is excited with a non-resonant 
frequency has an elliptic trajectory. The ratio between the 
semi-axes is given by the ratio between the frequency of the 
excitation and the resonance frequency. 

To test the applicability of the approximations leading to 
the analytical result in Eq. ( fT6] l we performed micromag- 
netic simulations for magnetic thin-film elements with differ- 
ent lengths, thicknesses, polarizations, and chiralities. The 
material parameters of permalloy are used, i.e., an exchange 
constant of A = 13 • 10~^^ J/m and a saturation magnetiza- 
tion of Ms = 8 • 10^ A/m. For the Gilbert damping we use a 
value of a = 0.01 which is in the regime as found by recent 
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FIG. 2: (Color Online) Dependence of the frequency uj and the 
damping constant F on the length I for various thicknesses t of the 
system. The symbols denote numerical results while the lines are fits 
with the analytical results. 



experiments 
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The degree of non-adiabaticity ^ is chosen 



to be equal to aJ^^ 

For the micromagnetic simulations we extended the imple- 
mentation of the Landau-Lifshitz-Gilbert equation in the Ob- 
ject Oriented Micro Magnetic Framework (OOMMF) by the 
additional current-dependent terms of Eq. ((TJiaiia^ The simu- 
lation cells are 2 nm in x- and y-direction which is well below 
the exchange length of permalloy. One cell of thickness t was 
used in z-direction. As in the analytical model we substitute 
the Oersted field by a homogeneous magnetic field. 

At first the four ground states with c ± 1 and p ± 1 are cal- 
culated for each / and t. The ground states are then excited 
by a short current pulse. The frequency ui and the damping 
constant F are obtained by fitting the subsequent free oscilla- 
tion with the first two terms in Eq. ( fT6l ). Results are presented 
in Fig.|2]and exhibit a good agreement between the analytical 
model and the micromagnetic simulationsiS 

For the driven oscillation we choose a magnetic film el- 
ement with length I ~ 200 nm and thickness t = 20 nm. 
This system size allows for reasonable computing time. The 
magnetization is excited with harmonic currents with a spin- 
polarized cuiTent density jP = 2.5 • 10^" A/m^ in ai-direction. 
The field excitation was performed with a harmonic field of 
H = 250 A/m in y-direction. The amplitudes and the phases 
of the oscillation in x- and y-direction of a vortex with positive 
polarization and chirality are depicted in Fig. [3] In the current- 
driven oscillation an excellent accordance between analyti- 
cal calculations and numerical simulations is found. In the 
field-driven case the amplitudes of the analytical solution are 
smaller than the amplitudes obtained from the micromagnetic 
simulations. These deviations are caused by the differences 
between the approximate magnetization depicted in Fig.[T]and 
the exact state. The phases between the maximum of the ex- 
citing magnetic field and the maximum deflection in x- and 
y-direction agree very well. Vortices with other polarization 
and chirality (not shown) yield the same accordancei^^ 
From Eq. ( fTSI l one can see that the current and field in- 
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FIG. 3: (Color Online) Amplitude of the (a) current-driven and (b) 
field-driven vortex oscillation in s-direction (solid red line, pluses) 
and j/-direction (dashed blue line, crosses) for a spin-polarized cur- 
rent density of jP = 2.5 ■ A/m^ and a field of H = 250 A/m. 
The insets show the phases between the maximum of the applied cur- 
rent or field and the core displacement in a;-direction (solid red line, 
pluses) and {/-direction (dashed blue line, crosses). The symbols de- 
note numerical results while the lines are derived from the analytical 
expression in Eq. l ll6t . 



duced forces on the vortex are of the same form. For experi- 
ments it is important to separate the Oersted-field and the spin- 
torque driven case. We describe the ratio between the field 
and current-induced forces on the vortex by tan = i^oe/^st, 
i.e., a mixing angle of C = and ( = ±7r/2 denote the 
fully spin-torque driven and the fully field-driven case, re- 
spectively. There are two possibilities to determine the ratio 
of both forces. On the one hand for non-resonant excitations 
the trajectory of the vortex core is elliptical as illustrated in 
Fig- SI According to Eq. ( fTSI ) the direction of the major axis 
of the ellipse is determined by C,. The amplitude of the vortex 
motion decreases very fast when the excitation frequency de- 
viates from resonance, i.e., for experimental observation very 
high current densities with frequencies close to resonance are 
needed. On the other hand the excitation mechanisms can be 
distinguished using the phase of the vortex deflectionjJi As in- 
dicated by the dots in Fig.|4]the position of the vortex at max- 
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FIG. 4: (Color Online) Analytically calculated phase between the 
maximum current or magnetic field and the a--deflection of the vortex 
core for a 200 nm x 200 nm x 20 nm permalloy square excited with 
a frequency of H = 4.8 GHz (above the resonance frequency of 
uj = 4.4 GHz). The inset shows a section of the sample with the 
simulated trajectories of the vortex core excited with i) (solid red 
line) a spin-polarized current density with an amplitude of jP = 
1.2 ■ 10^^ A/m^ and ii) (dashed blue line) a magnetic field with an 
amplitude of H = 1000 A/m. Points denote the position of the 
vortex at maximum current (i) and magnetic field (ii), respectively. 



imum current depends on C, which can be determined from 
Eq. ( fT6] l. The latter method is also applicable with excitations 
at resonance frequency. 

In conclusion we derived an analytical expression for the 
current- and field-driven trajectory of a vortex in thin-film el- 
ements. The analytical result is compared to micromagnetic 
simulations. The accordance between both approaches is very 
good. The analytical expression enables us to determine the 
ratio between spin torque and Oersted field driven motion. 
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